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Abstract 

We present a simple argument which determines the critical value of the anomaly 
coefficient in four dimensional conformal factor quantum gravity, at which a phase 
transition between a smooth and elongated phase should occur. The argument is 
based on the contribution of singular configurations ("spikes") which dominate the 
partition function in the infrared. The critical value is the analog of c = 1 in the 
theory of random surfaces, and the phase transition is similar to the Berezenskii- 
Kosterlitz-Thouless transition. The critical value we obtain is in agreement with the 
previous canonical analysis of physical states of the conformal factor and may explain 
why a smooth phase of quantum gravity has not yet been observed in simplicial 
simulations. We also rederive the scaling relations in the smooth phase in light of 
this determination of the critical coupling. 
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In recent work we have developed a framework which permits the study of the effects of 
gravitational fluctuations at large distances [|IJ-@- It is based on the quantization of a low 
energy effective action for the spin-0 or conformal part of the metric, which is determined 
by the trace anomaly. This action implies the existence of an infrared stable fixed point 
for quantum gravity in four dimensions with nontrivial scaling behavior |T| . On the other 
hand, quantum gravity can be reformulated as a statistical model on a random lattice and 
studied numerically ||. There is then the attractive possibility for testing the continuum 
predictions of an infrared fixed point by quantitative measurements in the lattice approach 

In a previous letter [[| we discussed the scaling relations for the partition function 
and observables in the conformal phase, based on the results of the infrared fixed point 
behavior found earlier. Since that time we have performed an extensive canonical analysis 
of the constraints of diffeomorphism invariance on the Einstein universe Rx S 3 , and found 
a discrete spectrum of physical states in the Fock space which survive the imposition of 
the constraints ||. Each of these physical states corresponds to an operator of scaling 
dimension or conformal weight 4 constructed from the conformal part of the metric and 
its derivatives. The existence of these states in a pure gravity theory without matter is 
in sharp contrast to the analogous quantization of the Liouville theory on the cylinder 
R x S 1 where there is only one state (the "vacuum") and only one operator (the identity 
or volume operator). In four dimensions, there is a tower of operators describing marginal 
deformations of the the theory away from its infrared fixed point. In the semi-classical 
limit, the two lowest of these become the Einstein-Hilbert and cosmological (or volume) 
terms. Each of these two operators comes with its own scaling behavior and anomalous 
dimensions, and each of them seems to imply a different value of the coupling where a phase 
transition to a highly nonclassical (or branched polymer) phase of quantum gravity could 
occur. Since the situation is more complicated than the D = 2 case, the scaling behavior 
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of observables should be reconsidered from the present vantage point of a more complete 
understanding of the physical states and the operators that create them, and the critical 
coupling determined by the point at which the scaling exponents first cease to be real. 

In two dimensions there is a simple argument for the critical value c = 1 of the em- 
bedding dimension in Polyakov's theory of random surfaces (or noncritical bosonic string 
theory) M. The argument is reminiscient of the Berezenskii-Kosterlitz-Thouless (BKT) 
argument for a phase transition in the 2D x-y model |J, in that it relies upon construct- 
ing a singular solution to the classical equations (called a "spike" in this context), which 
makes a contribution to the action that depends logarithmically on the infrared cut-off. 
Since the entropy of such configurations also grows logarithmically with the volume of the 
system, there is a critical value of the coupling at which the entropy always overwhelms 
the action and the partition function is dominated by a dense gas of such singular spikes. 
Conversely, if the coupling is adjusted in the opposite direction then such configurations 
are always suppressed in the infinite volume limit, and the geometry of the surface can 
be reasonably smooth. Since the coupling is also the coefficient of the trace anomaly in 
2D, this argument also tells us what the critical value of the matter central charge is, and 
why one should expect a phase transition from a smooth to a branched polymer phase at 
this critical coupling of c = 1. This phase transition has been verified in the dynamical 
triangulation approach to 2D random surfaces [f[D| . 

Since a branched polymer (or elongated) phase has been observed also in 4D simplicial 
quantum gravity || [7], it is quite natural to suspect that an analogous argument involving 
the "liberation of spikes" in the subcritical region of a coupling constant (Q 2 defined below) 
should apply again in this remarked in a recent paper |TT[ . In this note we develop 



this idea that the "half-wormhole" or "spike" configurations are the relevant ones for the 
phase transition at a certain critical value of the anomaly coefficient and we study the 
implications for 4D simplicial simulations and scaling relations in the smooth phase. 
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The Critical Coupling. In four dimensions, the analog of the Polyakov-Liouville theory 
of random surfaces is obtained by considering the effective action for the conformal factor 
of the metric induced by the 4D trace anomaly. This effective action takes the form [|TJ 

S E [a] = J d 4 x^g[aAta + \{G - |nJ2)a] , (1) 

where A 4 is the unique Weyl covariant fourth order differential operator on scalars, e 2a is 
the conformal factor of the metric (taken here with Euclidean signature), and Q 2 is the 
coefficient of the Gauss-Bonnet term G in the trace anomaly. It is normalized such that @ 

Q 2 = ^(N s + %N WF + Q2N V - 28) + Q 2 grav , (2) 

where Ns, Nwf, Ny are the number of free scalars, Weyl fermions and vector fields and 
Q 2 rav is the contribution of spin-2 gravitons, which has not yet been determined unambigu- 
ously. The —28 contribution is that of the a field itself which is the only known negative 
contribution to Q 2 . Thus, it is Q 2 which plays the role analogous to matter central charge 
and A 4 which plays the role of the kinetic operator □ in D — 2. In flat background 
coordinates A 4 = D 2 . 

Now we make the following observation. Because of the fourth order conformal differ- 
ential operator matched to the number of dimensions, the propagator of this operator is a 
logarithm, just as is that of □ in D = 2. This means that we have a situation analogous to 
that in 2D. Although the Mermin- Wagner theorem forbids a spontaneously broken phase 
with massless excitations, the 2D x-y model does exhibit a BKT phase transition at a 
certain critical temperature (i.e. coupling in the 2D Euclidean field theory), which is just 
the result of the logarithmic growth of the massless conformal invariant CP 1 propagator in 
D = 2. Because the conformal propagator Aj 1 in four dimensions has logarithmic growth 
at large distances, one should expect a BKT-like phase transition at a certain critical cou- 
pling Q 2 r in 4D for essentially the same reason as in the original x-y model or in the 2D 
theory of random surfaces. 
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Indeed, we can construct exactly the same spike solution in D = 4 as in D = 2, 

(Js(x) = q In \x — xq\ (3) 

with exactly the same interpretation. It is the solution of the classical equation following 
from (HO with a delta function singularity at X — Xq of strength q. In order to regulate 
the singular behavior of this configuration one may introduce an ultraviolet (UV) cut-off a, 
replacing \x — xo| 2 by 1 + \x — xq\ 2 /a 2 in the logarithm. This UV cut-off will be of order of 
the lattice spacing in the simplicial simulations. The infrared (IR) cut-off of the logarithm 
will be provided by the finite volume. 

Evaluating the action Se on this spike configuration in a large but finite spherical 
volume of radius L we find 

S E [a s } = 1 -Q 2 q 2 \n<^) . (4) 

On the other hand the number of ways of placing this configuration in the volume V is 
proportional to the volume, V ~ L A . Hence, the entropy grows like 41n(-^) and the free 
energy of such configurations behaves like 

*to]=&JV-4)h.(§) (5) 

for large L. It follows that for Q 2 > 8/q 2 these singular configurations will have positive 
free energy and be suppressed in the thermodynamic limit L — > oo, while if the inequality 
is reversed we must expect them to dominate the partition function. Consequently, there 
is a critical value of the coupling Q 2 at which we expect to see a phase transition from 
a phase with many sharp elongated spikes to one where such singular configurations are 
suppressed. The analogy to the BKT argument involving vortices in the x-y phase field is 
obvious. 

In two dimensions, Cates argued that the value q — —1 (in the present notation) is 
special because this is the strength at which the singularity is first strong enough to give a 
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divergent contribution to the classical volume (area in D = 2) as the cut-off a is removed 
||. Hence, these are the configurations with the lowest free energy that can dominate 
the partition function in the continuum limit L ^> a. Since the volume or cosmological 
operator is a marginal deformation of the free Liouville theory, this argument for q = —1 
can be turned into a renormalization group analysis in the IR as well, by taking into account 
the gravitational "dressing" of the volume operator due to the loop corrections of the free 
□ propagator near its Gaussian fixed point. Then we find that q = — 1 is precisely the 
condition that justifies the neglect of the cosmological term and the use of the spike solution 
to the free Liouville theory in the infinite volume limit L — > oo. Let us give this argument 
why q = — 1 is the relevant spike solution in four dimensions as well. 

In four dimensions, there are two terms with fewer derivatives that can be added to the 
free action (|l|). They are the volume or cosmological term, 

\S [a] = A J d 4 x^ = A J d 4 xe 4a -> A J d 4 x . (6) 

and the Einstein-Hilbert term, 

□ (da) 2 



- -— f d 4 x 
Ak J 

The integrands in these expressions have engineering dimensions and 2 respectively, and 
so they must be multiplied by powers of e a to give a scalar density with total conformal 
weight 4. In the last forms of (||) and (0) we have allowed for the possibility that these 
terms are gravitationally dressed and that the classical scaling codimensions (flo) c i = 4 
and (/?2)c/ — 2 are modified at the quantum level. For the Einstein term involving two 
derivatives of the a field there is also the possibility of renormalization group mixing with 
operators of lower dimension, which is represented by the f{Q 2 ) term in eq. (0). Indeed, a 
covariant computation of the renormalization of these operators due to the loop effects of 



TT^M = -7T / d 4 x^R = - / d 4 xe 2a 



72tt' 



K 



f(Q 2 



(7) 
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the free D 2 propagator shows that [3q and satisfy the equations, 



A = 4 + ^, (8) 



and 



2Q 2 ~ 

while the operator mixing is given by 



& = 2 + £%, (9) 



2, .. " 2 



HQ, Q2 



Aa 2 6a 4 
1 + ^rr + 



(10) 



Q 2 Q 4 

and a = ^2/2 in the present notation The classical values for the codimensions (3q and 
fa are recovered only in the limit Q 2 — » 00, where quantum fluctuations of the conformal 
factor are suppressed. Notice that this classical limit is the analog of the opposite limit of 
the central charge (i.e. c — > —00) in the 2D case, since Q 2 is positive for free conformal 
matter fields. 

Corroboration of the anomalous scaling dimensions ([8]) and @ comes from the canonical 
quantization of the a theory on the Einstein space Rx S 3 where an infinite tower of discrete 
diffeomorphic invariant states, labelled by the integers, was obtained ||. In the semi- 
classical limit Q 2 — > 00, these states are created by operators which are volume integrals 
of integer powers of the Ricci scalar, i.e., 



J d 4 x v^iT ~ J Ae fe " ff [(da) 2n + •••], (11) 



where 



f3 2n = 4-2n+§%. (12) 



2Q 2 

The ellipsis in (|TI|) refers to the operators with lower numbers of derivatives which mix with 
(da) 2n under the renormalization group at finite Q 2 , analogous to (0) and ( |T0"D for n = 1. 
The exact forms of this operator mixing is not determined by the canonical analysis, but the 
values of all the /?2n are fixed by the Hamiltonian state condition (the analog of the Lq — 1 



-8- 



condition on the physical states in D = 2). The first two of the discrete physical states are 
created by operators of the form (|6|) and (0) and the values of the scaling exponents are 
identical to (§) and ([|) obtained in the covariant Euclidean approach. 



Now the volume term evaluated on the spike configuration behaves like 



So[as] 



(4 + A,g) 



c - - 
a 



2a 4 Q 2 



L 

c - ( - 

a 



f4 

2Q 1 



(13) 



for L>a and q = —1, which is the only value for which the integral is IR convergent for 
any positive value of Q 2 (for which Re (3q > 4). In this expression cq is a positive constant 
of order unity whose value depends on the precise way the UV cut-off a is introduced. In 
a similar manner, the Einstein term evaluated on the spike configuration behaves like 



S2 Ws] 



q 2 a 2 



(2 + /3 2 q) 



2+p 2 q 



2a 2 Q 2 

' 0i 



L 

'--2 - [ - 
a 



4 

2Q 2 



(14) 



for L 3> a and q = —1, where C2 is another constant (dependent on a 1 Ik). Inspection of 
the previous form shows that q = — 1 is again the only value of q for which the Einstein 
action evaluated on the spike configuration is convergent as L — ► 00, for any positive value 
of Q 2 (for which Re (3 2 > 2). 

Hence, knowledge of the scaling exponents of the Einstein and volume operators allows 
us to compute their values on the singular spike configurations and show that they are 
subdominant to Se in the infinite volume continuum limit for all sufficiently large and 
positive Q 2 , if and only if q — — 1. With q — — 1 selected in this way, we obtain from our 
previous evaluation of the free energy of the spike (§) the critical value, 



Q 



(15) 



For Q 2 > 8 the solution of the quadratic relations (H) and (||), viz. 



(16) 
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and 

h = Q2 { l -f-?) (17) 

are indeed both real and Re (3q > 4, Re (3% > 2 for 8 < Q 2 < oo, so that the previous 
argument for the irrelevance of the Einstein and cosmological terms in the infinite volume 
limit are justified a posteriori. We have chosen the minus sign for the square roots in (|Tjj) 
and (|T7|) in order that /3 2n — > (/?2n)d — 4 — 2n in the semi-classical limit Q 2 — >■ oo. 

When Q 2 < 8, then one of these exponents (namely /?o) becomes complex and we should 
expect qualitatively different behavior of the theory. This is an independent indication of 
the critical value Q 2 r = 8. From the simple BKT-like argument, we expect that Q 2 < 8 
corresponds to a phase which is dominated by a dense gas of q = — 1 spikes. The metric of 
the singular q = — 1 configuration (with the origin, x = 0) is 



ds 2 = e 2as ds 2 = ^ + J 2 2dQ2) = (rflnr) 2 + dQ 2 



which is just that of R x S* 3 , if we start with the flat R 4 metric background, ds 2 . Since 
the S 3 has arbitrary radius and the R is the entire real line, this geometry is arbitrarily 
elongated and thin. Since any region is locally like that of R 4 , a "gas" of such configurations 
means that these "spikes" (or perhaps more appropriately, "tubes" or "punctures") will 
look like a large number of branches with a hierarchy of thinner and thinner sub-branches. 
A random geometry with many such spiky extrusions looks very much like the elongated 
phase described in ref. [^]. Conversely, if Q 2 > 8 then the scaling dimensions (|1^) and (|T7|) 



are real, the spikes are suppressed and one would expect to find that the partition function is 
dominated by much smoother configurations, which at very large Q 2 are well-approximated 
by semi-classical metrics. 

Thus, Q 2 r = 8 behaves in many respects like the c = 1 case in 2D gravity, with the 
important difference that the addition of conformal matter fields brings us into the smooth 
phase Q 2 > 8 rather than into the elongated or branched polymer phase c > 1 (see eq. (0)). 
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Now in pure simplicial gravity, i.e. without the introduction of matter, only an elongated 
phase and a crumpled phase (where the geometries collapse upon themselves) have been 
found. There is a transition between them as the lattice coupling corresponding to the 
Einstein term is varied, although apparently there is as yet no universal agreement as to 
the order of this transition. From the point of view of the considerations presented here the 
numerical results are consistent with the hypothesis that pure gravity has Q 2 = Q 2 grav ~ 8, 
and only a spiky elongated phase in the continuum limit. We should emphasize in this 
connection that the anomaly generated action Se is positive definite for Q 2 > so that there 
is no conformal factor problem, nor any need for a conformal rotation as in the Einstein 
theory. Hence, the statistical continuum theory of random four-geometries described by 
(0) is well-defined for < Q 2 < 8 even if these geometries turn out to be very jagged. 

The main theoretical difficulty in determining Q 2 rav is that the Einstein theory is neither 
conformally invariant nor free, so that a method for evaluating the strong infrared effects of 
spin-2 gravitons in the continuum must be found that is insensitive to ultraviolet physics. 
In ref. || we performed a strictly perturbative computation which gives the value Q 2 rav = 
1411/180 ~ 7.9 for the graviton contribution. Since the method used is based on the heat 
kernel expansion, there is the possibility of mixing up ultraviolet with infrared effects in 
this evaluation, and we cannot regard it as definitive. However, a computation using the 
totally different conformally invariant Weyl tensor-squared action leads to a similar value 
Q 2 grav = 8.7. Hence it is likely that the correct infrared graviton contribution to Q 2 is in the 
neighborhood of 8, and that most of this large value is due to the spin of the field, which 
is larger than any other individual field's contribution of lower spin. If the value of Q 2 in 
the pure gravity theory turns out to be less than eight, then one would not expect to find 
a second order phase transition to the continuum limit in the simplicial simulations with 
geometries that are smooth. Instead the continuum limit of the random geometries would 
be an elongated phase, filled with many spiky extrusions. 
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If this is indeed the correct interpretation of the numerical results to date, then the 
introduction of some number of free conformal matter fields to the simplicial simulations 
would push the value of Q 2 above 8 and lead to a continuum limit exhibiting a smooth 
phase, perhaps by softening a weakly first order transition between the elongated and 
crumpled phases into a second order transition with the new phase appearing in the phase 
diagram extended into the Q 2 > Q 2 grav direction. Taking the 1411/180 number at face value 
and using the known contributions to Q 2 for lower spin fields in eq. @, one finds that it 
would take the introduction of 57 conformal scalars, but only one photon to induce this 
transition to the smooth phase. 

Scaling relations. Assuming that the smooth phase with Q 2 > 8 exists, the exponents 
make definite predictions for the scaling of the fixed volume partition function, 

Z(k, A) = J [Da] exp [-S E [a] - ±-S 2 [a] - ASqM) • (19) 

The covariant continuum measure [Da] has been discussed in several previous articles 
0, 0. By inserting a delta function of physical four- volume, we obtain the fixed volume 
partition function 

Z(k; V) = j [Da] exp (-S E \a] - ^S 2 [a]^j 5 (S [a] - V) . (20) 

If the effective action appearing here, namely the sum of the Einstein and anomaly 
induced terms, is written in an arbitrary curved background and one studies the effect on 
the finite volume partition function of a shift in <r, 

a^a + uj (21) 

then one finds that 

Z{k- V) = exp \-u (A) + Q 2 X E )] Z(Ke'^-e-^V) , (22) 
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where the Euler number x E = 2 for fixed S* 4 topology. Since a has been integrated out, 
the resulting Z(k; V) must be independent of the shift (0), which is only consistent with 
(E2|) if the fixed volume partition function can be expressed in the form, 



02 . 



Z(k; V) = V~ X ~ 2 ^ Z(kV~^) , (23) 



for some function Z of a single argument. This scaling relation differs from that obtained 
in eq. (10) in our previous letter since now we have allowed the volume and Einstein 
terms (§) and (|7|) to have independent scaling exponents in agreement with our subsequent 
canonical analysis |j) and ref . 1 13 . 



By using the relation (|I6|) we now obtain the susceptibility exponent 



7(Q 2 ) = 2- 2 g = -2 V 1 (24) 
Po 1 - Jl - 



Y2 



instead of eq. (23) of ref. ||. Notice that if Q 2 is close to 8 then the susceptibility exponent 
is close to zero, which is what has been found in the simulations to date. As shown in 
Fig. 1, 7 is negative for Q 2 > 8 and approaches zero from below as Q 2 —>■ 8. 

The form (|23|) also implies that we should require kV o to remain finite in the scaling 



limit, or equivalently, 



with 



k~V s as V oo , (25) 

02 _ 1 " i 1 



1 

Q- 



S=^ = \ . (26) 



We remark that this relation differs from eq. (22) of ref. 0] because of the addition 
of the second independent cosmological operator which we did not consider in that work. 
The exponent 5 is plotted as a function of Q 2 in Fig. 2. For Q 2 near 8, 5 ~ 0.3, which 
differs from the value 0.47 ± 0.03 reported in the simulations of ref. fTj. However, if our 
hypothesis is correct, namely that Q 2 for pure gravity is less than 8, then the smooth phase 



Fig. 1 




Fig. 2 




Figures 1 and 2. The susceptibility and scaling exponents 7 and 5 as functions of the 
anomaly coefficient Q 2 . 



where eq. (Eq) applies has not yet been produced in the simulations. Of course, when the 
numerical situation has been clarified and an unambiguous continuum limit with a smooth 



phase has been demonstrated, then there is only one free parameter and both (p4) and 
must be consistent with the same value for Q 2 , if the conformal fixed point relations 
predicted by the continuum theory are correct. 

Finally, one can consider any operator composed of either metric or matter fields with 
scaling dimension A, in the absence of dressing. By inserting this operator in the fixed 
volume partition function 



(Oa) v = J [Va] O k exp (-S E [a) - ±-S 2 [a\) 8 (S [tr] - V) 

and repeating the shift ([H]) and scaling argument for this quantity, we find 



(27) 



(2* 



where the full scaling dimension A is related to the "classical" dimension A by 
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Here is the relevant codimension of the gravitational dressing determined by 

/?A = 4-A + ||. (30) 

At large Q 2 one has: 

A - A = s| A < 4 - A > = ^ A < 4 - A > + '" < 31 > 

For the volume and Einstein operators for which A is and 2, respectively, we recover the 
previous relations (||) and @. 

If these scaling relations can be verified in the smooth phase, and the infrared conformal 
fixed point predicted by the action (|I]) confirmed, then the lattice could provide a nonper- 
turbative method for measuring the contribution Q 2 grav . In addition to being interesting 
in its own right by exploring a nontrivial fixed point of 4D quantum gravity, the scaling 



relations could also find applications in cosmology [JT4 . 
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